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Abstract 

fl} A theory of the t-J model in the presence of an external electromagnetic field is presented. The 1/N 

P^ expansion for this model in the slave-fermion representation is developed, and it is shown that all the 

general properties of the theory are satisfied in every order of the 1/N expansion separately. A convenient 
procedure of the gauge fixing is suggested. It is shown that superconductivity may exist in this model even 
in the absence of the electron pairing. It is argued that experimental predictions of the theory might agree 
with the observed properties of the copper oxides. 
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^ I. Introduction 



Soon after the discovery of superconductivity in the rare earth-based copper oxides [1], Anderson [2] suggested 
*^3 that the physics of these materials may be described by a large-?/ Hubbard model and the related t-J model. 

The Hamiltonian of the t-J model operates in a constrained subspace which excludes double occupancy. To 
cope with the constraint, the electron operators are usually represented as 



4r = 4r«r- (1.1) 

The operator a r here represents an empty site, while the operator b aV represents the electron's spin as S r = 
\h aY u a jjb jjy ■ In terms of these operators, the constrained subspace is selected by the holonomic condition 

bl r b ar + ala r = 1, (1.2) 

which has to be imposed independently at every site r. Representation (1.1) simply expresses the fact that 
the states of the electron system with forbidden double occupancy may be put into one-to-one correspondence 
with the states of the b-a system, satisfying condition (1.2). The corresponding map is given by Eq. (1.1). 
To preserve the fermionic statistics of electrons, we must assign bosonic statistics to one of the particles a or 
b and fermionic statistics to the other. Either possibility is admissible, and we choose to assign the fermionic 
statistics to a; the particles b then are usually called Schwinger bosons (SB). 

Representation (1.1) corresponds to the separation of spin and density degrees of freedom [3]. Their 
excitations - spinons and holons - rather then electrons, have been claimed to be the real quasiparticles in 
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the copper oxides. If this is true, then what might be the nature of the superconductivity in the system of 
interacting spinons and holons? 

The phenomenon of superconductivity may be described as an emergence of the photon's mass due to the 
Anderson-Higgs mechanism. The Higgs field, which screens the photon, is necessarily bosonic; to create such 
a field, electrons in the conventional BCS theory have to form a bound state - a Cooper pair. On the other 
hand, in the spin-density separation approach to the copper oxides, we have the boson field b a as one of the 
fundamental fields of the theory. Then the natural question to ask is whether this field can screen an external 
electromagnetic field (EMF), leading to superconductivity without electron pairing? 

It is clear that the only way to answer this question is to develop a reliable approach to the t -J model in the 
presence of an EMF, and to check whether an external magnetic field in this model can create a nonvanishing 
electric current in the phase having unpaired electrons. This program is realized in the present paper, and we 
demonstrate that the question raised above has a positive answer. We note that the spin-density separation 
cannot be considered a trivial consequence of representation (1.1). Indeed, this representation naturally 
introduces to the theory invariance with respect to the local gauge transformation a r (& CTr ) — > e 8a i* ar (6 (7r ). 
Thus the theory becomes a lattice gauge theory [4], where the presence of the fields b a and a in the initial 
Lagrangian does not automatically guarantee their existence as real particles. Indeed, if the gauge field is 
in the confining phase, it binds spinons and holons together to form a physical electron as a bound state. 
The electron then carries both spin and density degrees of freedom, so that in effect there is no spin-density 
separation. 

The real spin-density separation occurs as a result of spontaneous gauge symmetry breaking with the 
transition to the phase with nonzero average gauge field and unconfined slaves b a and a (deconfinement 
transition). That kind of transition was discovered on the mean-field (MF) level in the very first works on 
lattice gauge theory [5, 6]. Soon, however, Elitzur [7] proved that in the lattice gauge theory the average value 
of the gauge field must vanish. The reconciliation, achieved in numerous later works, may be described as 
follows. Initially, one does not need to fix the gauge in lattice gauge theory [5], because the gauge group there 
is compact. In such a theory without gauge fixing and with manifest gauge invariance, the expectation values 
of the gauge-noninvariant quantities do indeed vanish due to the Elitzur theorem. On the other hand, one 
can fix the gauge, using the Faddeev-Popov trick, and factor out the trivial volume of the gauge group. This 
procedure (which is also necessary for the development of the diagram technique) will not affect the physical, 
gauge invariant results. Fixing the gauge, however, spoils the gauge invariance of the theory. The Elitzur 
theorem is not applicable to it anymore, and the gauge-noninvariant quantities can obtain now nonvanishing 
expectation values. Therefore, on one hand, the lattice gauge theory with the fixed gauge gives correct values 
of the gauge-invariant quantities, and on the other hand, allows one to give a convenient description of the 
different phases of the theory. Indeed, here we can use for that description the expectation values of the gauge- 
noninvariant quantities, such as the gauge field itself. This is more convenient than the use of gauge-invariant 
(and more complicated) objects, such as the Wilson loop [4, 5], in the theories without gauge fixing. 

In this paper we will follow this route in the study of spontaneous gauge symmetry breaking and spin-density 
separation. The outline of the paper is as following: we start in the next section with the derivation of the 
diagram technique for the t-J model. We develop a systematic 1/N expansion, and demonstrate throughout 
the paper that all general properties of the theory are satisfied in each order of this expansion separately. The 
zeroth order, corresponding to the N — > oo limit, yields the mean-field theory (MFT). We work consistently 
on the lattice without the transition to the continuous limit and suggest a convenient procedure of the gauge 
fixing. In the following sections we discuss the general properties of the theory: Sec. Ill contains the discussion 
of collective modes, which exist in the phases with a broken gauge symmetry, while Sec. IV shows the relation 
between the diagram technique and the equations of motion, following from the Lagrangian. The interaction 
with an EMF is introduced in Sec. V. In Sec. VI we demonstrate that this interaction satisfies the Ward 
identities. The superconducting phase is considered in Sec. VII, and finally, in Sec. VIII, we discuss the 
possibility of describing the experimentally observed properties of the copper oxides in the framework of our 
approach. Since this paper is devoted mainly to the formulation and the basic consideration of our approach 
to the t-J model, the explicit comparison between our theory and experiment will not be reported here. 



II. 1/N expansion and diagram technique 

II. 1 Basic definitions 

Our starting point is the t-J model, which can be represented by the Hamiltonian [8] 



r v 



1 (~<\ r< 



(2.1) 



where ijj' aV (Vvr) are the electron creation (annihilation) operators with spin a at site r (here r is a vector with 
integer- valued components), r\ runs from 1 to the dimension of the space d, r\ is a unit vector in a corresponding 
direction, and C v r = V'trV^r+ij — VM.rV'tr+jj- The second term in (2.1) is just the superexchange interaction 
J v (Sr • S r _|_,j — jTi r nr+fj) , where S r is a spin operator and n r is a number operator at site r. 

In (2.1) we have taken into account that high-temperature superconductors are all strongly anisotropic 
materials with direction-dependent parameters. To avoid the countless repetitions of the sign of summation 
over T), we will utilize an agreement (the generalization of the usual Einstein rule) that this sign will not be 
written explicitly in formulas where the index r\ is used two or more times independently of its position. For 
example, summation over r\ will be assumed in such expressions as t~ 1 e tkr > cos k v . The geometrical meaning of 
this rule becomes clear in a long-wavelength limit, where all expressions of that type turn out to be scalars, 
created in the usual way from vectors and diagonalized tensors. This rule will be not applied in cases when 
the expression with repeating ry-indices stands on one side of an equality, if the other side of the equality also 
has the index rj, but used only once. 

As was pointed out in the Introduction, in order to impose the single-occupancy constraint, the electron 
operators are usually represented as a product (1.1). We will see however, that it is more convenient, especially 
in working with diagrams, to use a different, although equivalent, representation 

4r = blrfl (2-2) 

where /jl = a r . The constraint now has the form 

btrKr = flfr- (2.3) 

Due to the fermionic nature of the operators / r , the constraint (2.3) automatically restricts the possible number 
of SB and, of course, electrons at site r to or 1. 

The representation similar to (2.2) with the constraint (2.3) works equally well for the large-?/ Hubbard 
model with more then one electron per site (i.e., for electron doping). There, the constrained subspace has 
one or two electrons on every site; the states with empty sites are excluded. We can then consider a doubly 
occupied site as corresponding to the absence of b and / fields, and a site with one electron as corresponding 
to the presence of one / and one b particle. An empty site would correspond here to the presence of two 
/ particles at one site, which is again impossible due to their fermionic nature. In each case, therefore, the 
presence of an / particle means the presence of one and only one electron, so we will use for them the name 
monon. 

The physical interpretation of the representation (2.2) is that in order to satisfy the condition of no doubly 
occupation, we can choose the basis many-electron wave functions of the constrained space in a form 

* CTl ,...,^ (ri, . . . , rjy) = & (n, . . . , v M )^ b au >aM (n, . . . , rjy), (2.4) 

where ^ (^ b ) is an antisymmetric (symmetric) function of its arguments. We can then regard excitations b 
and /, which correspond to wave functions 'J 6 and ^ , as the original elementary excitations in our theory. 
Their occupation numbers must satisfy the constraint (2.3), and their relation to electrons is given by (2.2). By 
(2.4), this relation is local in a sense, that both particles b and / occupy the same sites as electrons. This is in 
contrast to the nonlocal relation, corresponding to the usual representation (1.1), where spinons b occupy the 



same sites as electrons, while holons a occupy the remaining sites. It is interesting to note that the slave-boson 
formulation, which uses the same representation (1.1), but assigns fermionic statistics to operators b aV and 
bosonic statistics to operators a r , cannot be transformed to the monon language. 
The substitution of (2.2) into (2.1) gives on the constrained space 



^ = E 



-t v [B nY F nY + B vJ .F vJ .) —A vr A v r 



(2.5) 



where B v r = bij.brjr+rj, F v r = frfr+rj and A v r = bfrbir+rj — birbfr+rj- For the spinon's fields we will use 



Nambu's notation: 

/ u-fr 



f ) , bt= (&} r , b ir ) (2.6) 



and Pauli matrices Tj , j = 1, 2, 3 and t± = T\ ± ir^, besides, we denote ir± = 1 ± T3. For the operators A v r, 
B v r, F vr and their conjugates we will also use the general notation M avr , a = 1, . . . , 6 with the identification 

M av r = [A vr , A v r, B vr , B vr , F vr , F vr \ . (2.7) 



We have then 



M li2r ,r = ±- f4 r ± & r+ij - b\ + -T±b Y 



M-s, MY = - \blTT T b r +^ + b[ + -TT ± b x 
Besides operators (2.7), we also introduce the operator 



(2.8) 



M 7V = i (b\b v - fl f v y (2.9) 

We will denote all M-operators as M avr , keeping in mind that M^ does not have the index r\. For the 
frequently used sums over a we will utilize the following notations: ^2 a = Ea=i> E« = X^aeA^ ' anc ^ 
E« b>/ =E ae ^ / ,where^ 6 = {l,...,4},^ / ={5,6},A ? -6 = {l,...,4,7},andA ? > = {5,6',7}. 
The constraint has the form 

bUr = fifr, (2-10) 

or M^r = 0. In all formulas which include the product of SB operators in Nambu's notation, the normal 

ordering of initial operators & CTr , &^ r is assumed. 

A convenient way to work with Hamiltonian (2.5) and the constraint is the 1/N expansion. We generalize 
our model to N > 1 colors by setting t v — > t v /N , J v — > J v /N , &r(/r) —* &r (/r )> a = 1, • • • , A. The product 
of boson or fermion fields within each composite operator M avr is now understood as a scalar product with 

respect to the color index; the constraint (2.10) is generalized to &r &r = /r /r • I n ^ ne following, we will not 
write the color indices explicitly, since it will be always clear what kind of scalar product with respect to them 
is used. 

To impose the constraint, we add to the Hamiltonian the new term H e = — (2Ae) _1 J^ r My r . As was 
demonstrated in [9], keeping e finite helps to control infrared divergences, which must cancel in each order of 
the 1/A expansion separately. The limit e — > must be taken at the end of the calculations. We now introduce 
the 7x7 matrix h v with nonzero matrix elements h n Yi = h V 2i = </jj/2, h v s5 = h V 5s = h V 46 = /*jj64 = ^jj, 
h v 77 = 1/e and write the partition function as an imaginary-time path integral: 



Z = V(b*bf*f)exp - dxlb* x T 3 b x + f*f x 

JT7 2_^ ^arixhriafjMfjj^x — fj,n x 



2A 

a/3 



Here the chemical potential jjl fixes the electron density and we have used the notations n x = (l/2)(b x b x + f* f. 



r/ 3 



x = (t, r), J dx = L drJ2 r . We now introduce the real auxiliary field (AF) <f>? x and the set of complex AF 
b artx , a = 1, . . . , 6, which are divided into three pairs of conjugated fields according to 



'2i;a; 



'li)ii Y'srix 



a 3i)j;) <P6!j:c 



"Atix- 



(2.12) 



We will 



for all AF the notation 



> ar]x 7 



with the understanding that (f>r x does not have the index r\. Then 



the interaction terms in (2.11) may be decoupled, yielding 



Z 



V(b*bf*f^)exp 



dx( b x T 3 i> x + f*f x + - ^2 <?Wn^° c 

J ^ a/3 



afiVPnx 



/ j YOiTjx^arjx t^^x 



(2.13) 



(o) 



where the bare polarization matrix of AF n}j is given by 



n 



(0) 
rja(3 



N(h- r ) 



V J a/3 



for h val 3 = 

N 
■ for h val3 jt 



(2.14) 



Equation (2.13) is invariant with respect to the local gauge transformation, induced by the representation 
(2.2) and the constraint (2.3). To make this gauge invariance explicit, we use the following parametrization of 

,6: 

<t>a-qx = \<t>a-qx\ exp { U [g a ar, x + g a {a r]x + jaQrix)]} , ( 2 -15) 



AF <j> ar , x , « = 1 



where / is a lattice constant and g a , g a , j a , a = 1, . . . , 7 are given by 

g a = (-1, 1, 0, 0, 0, 0, 0), 



(0, 0, 1, 
(0, 0, 1, 



-1, -1, 1, 0), 
-1, 1, -1, 0). 



(2.16) 



We have then 



Y'i'Yli9a^(t> aV x, 

a 
JTj^2galn<l>a V x- 



41 

a 

Now in the limit e — > Eq. (2.13) is invariant with respect to the local gauge transformation 

h v e iXx b f v p-^'x* f 

u ax 7 c u ax i J x 7 c J x i 

®0x — > ®0x + 9qXx, 

@"nx / @"nx \ 0„ \x i @"nx ' @"nx \ GL Xx • 



(2.17) 
(2.18) 

(2.19) 
(2.20) 
(2.21) 



Here ao x = <f>r x , do = d/(dr), d^ is a "right" lattice derivative, d^Xx = l~ 1 {Xx+rj — Xx), where we define 
x±f) = (r, r±fj) and d?' is a "staggered" right lattice derivative, defined as d?' = £,rd^£,r- Here £ r = exp(i7?r), 
where 7? is a vector which has all components equal to tt; obviously, d^Xx = ~l~ 1 {Xx+rj + Xx)- We see, that 
the gauge invariance in this problem has a nonstandard form, being characterized by the presence of two gauge 
fields: "staggered" d vx and "uniform" a vx . 



and self-consistency equations 



II. 2 Diagram technique 

II. 2.1 Real space 

Due to the gauge invariance, Eq. (2.13) contains redundant integration, and in order to develop the diagram 
technique we must fix the gauge. We will do this by adding to the Lagrangian in (2.13) the gauge-fixing 
term (GFT) C g {x). The explicit form of the GFT will be presented later. After fixing the gauge the original 
local symmetry is lost, so the remaining symmetries are global and may be spontaneously broken. Due to 
this spontaneous symmetry breaking, the spinon fields b x as well as AF (f) avx and composite fields M arjx , 
a = 1, . . .,6 may obtain nonzero expectation values. As was discussed in the Introduction, although these 
fields are not gauge-invariant, with a fixed gauge this is possible and does not contradict Elitzur's theorem 
[7]. We also assume the possibility of having nonzero expectation values for the field (f>r x and the composite 
operator Mj x (see the corresponding discussion below, after Eq. (2.27)). Consequently, we make in (2.13) the 
shifts b x — y v r + b x and (f) avx — > <j) art + 4> avx , where v r and <f> ari are the expectation values of the corresponding 
fields while b x and (f) avx are now their fluctuating parts. Here parameter v r is time and color independent 
and for the general discussion in this paper we only consider a simple case of space and time independent 
parameters <j) ari . The choice of v r and <f> ari will be specified later. It will be clear then, that these parameters 
are determined on the MF level and do not depend on the gauge choice. Therefore, it is legitimate to choose a 
gauge condition, which itself depends on these parameters. In other words, we first determine the phase of the 
model, and then choose a gauge, which is appropriate in this phase. At this moment we will leave the GFT in 
the representation (2.22) below in an unspecified form. 

We now arrive at the final path-integral representation for the partition function 



Z 
£{x 

Cmf{x 



C 



(0) ( 
AF\ 



t^int \% 



V(b*bf*f^)exp 



C{x)dx 



,(0) 



£mf{x) + jCn n (x) + C A ' F (x) + C g (x) + £i n t(x) + const , 

K T 3bx + f x fx - /2^ a, l Ma, l x ~ ^ Ux) 
a 

-^2<t>a^M a ^ x - Hb(b* x v r + v\.b x ), 

-Td, n (0) 

9 / .fayx'-'-na. 



(2.22) 



qafj'rP'nx, 



aj3 



■ ^2 ^a VX M' a 



ar]x 7 



where 



M n 



M' 

arjx 



M n 



M n 



M, 



at] 



l^£ni 



yv-^-\-b x 



yv-£-\-b x 



M n 



NM n 



r(0) , 



M n 



b x ^-Vy 



(2.23) 



and we denote /^j = ±i<f>7 + A«/2. Cu n (x) includes terms, that are linear in spinon fields. 

The representation (2.22) generates the diagram technique in our theory, with Cmf regarded as a free 
Lagrangian of b and / fields, which determines their bare propagators. On the diagrams, we will represent 
the propagators of fields b and / by solid lines, and the propagators of AF (which will be determined later) 

by wavy lines. The factors v r , Vj. will be represented by short zigzag tails. We note, that if we express 
Cmf as 



' ar] 



-xri — 2_^ h^cfiM/ 



Pv 



(2.24) 



then Cmf will correspond to the MF decomposition of the Hamiltonian (2.5). 

In the following, we will use for the arbitrary functional of the fields T the notation 



{T) = Z~ l V(b*bf*f<f>)Fe 



■ I C(x)dx 



(2.25) 



The notation (T)mf will have the same meaning, with the changes C — > Cmf and Z — > Zmf, where Zmf is 
the partition function of the MFT with the Lagrangian Cmf- 

To determine the unknown parameters //, v r and <f> ari (or M arj ) we will use the following conditions. First, 
if S is a density of holes, then we must have 



(n x ) = N(l-S). 



(2.26) 



Second, we require that the part of the action linear in b x , b* x , i.e., Cu n {x), must vanish; it will become clear 
below that this condition yields limjv->oo(&a;) = 0. We will return to the explicit form of this condition later. 
Finally, we require the vanishing of the "bare tadpole" subdiagrams of Fig. 1(a), so that the diagrams with 
these subdiagrams must be excluded. This condition is (M' x )mf = 0, or 



(M< 



ar]x 



IMF 



NM, 



b x — >-vj*-\-bx 



QL'q ■ 



(2.27) 



Thus the MF expectation values of operators M ar]x per color must be equal to the values M ar) , introduced in 



(a) 




where 




NM, 



a i] 



(b) 




FIG. 1. Forbidden (sub)diagrams. (a) "Bare tadpoles." The shaded circle represents {M ar)x )MF = {M ar)x )MF — N M ar) . 
Shown is the case vr = 0; when vr ^ and a 6 Mb, the obvious term with two zigzag tails must be added to the 
shaded circle, (b) RPA bubble, inserted to the AF line. 



Eq. (2.23). For a = 7 we note, that with our softened way of introducing the constraint, the operator M^ x is 
allowed to have a small expectation value, controlled by parameter e. In order to obtain the finite results in 
the limit e — > 0, we take ^7 to be independent of e, so that M7 = £^7 vanishes as e for f->0. The value of <f>^ 
itself is then determined by Eq. (2.27) with a = 7, which in the limit f->0 gives 



\v r \ 



(Kb 



x IMF 



\Jx Jx 1 



MF , 



(2.28) 



where we used the equality (b x )MF = 0. We note that the parameters M arj are determined on the MF level 
by equations (2.26) - (2.28); the exact expectation values (M avx ) might be, of course, different from the 
corresponding MF expectation values, used in these equations. 
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II. 2. 2 Momentum space 

We will now derive explicit expressions for the bare Green's functions of spinons and monons and their in- 
teraction vertices with AF. We denote k = (w,k), kx = uit — kr and J2 k = (j3V)~ l ^ w ^ = T^2 W j„ ,„ , d , 
where w is a Matsubara frequency, k is a lattice momentum, V is the number of lattice sites, and £1 refers to 
the first Brillouin zone. The Fourier transformation is defined according to 



ipk= dxip x e ikx , i> x =^2ip. 



fee 



-ikx 



(2.29) 



for an arbitrary value ip x , while the value ipl is transformed according to the conjugate definition. For the 
Nambu spinors (2.6) this gives 



°i-k J 



The Fourier components of the operators M are 



M, 



ar/q 



e' 2 



^a^k + i 



eM b 



& ^r 4 



M. 



79 



fe 

J2 ( h k-i h k+i-r k _jk 



lk+i 



(2.30) 



(2.31) 



where 



m 5,6„k 

For the MF action in Eq. (2.22) we have 



m i,2i,k = ±*V±sin^, 
1 
2 






m 3,4„k = 7^T e ""'+ 7r ±e 



IK n 



>MF 



E 



5 fe -«wr 3 - //{, 



X^ m a 



+ /fc ~ iuJ ~ VS ~ J2 <^ m a„k fk 



The bare Green's function of spinons, G x = (b x b* Q )MF , has the Fourier transformation 

_ 1 1 ( (hk h \k)MF (b^kbl-k)MF \ 



From (2.33) we have for it 



Gk 



-IUJT 3 - fl b 



X^"W 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



The bare Green's function of monons g x = (f x fo)MF nas the Fourier transformation 



1 



9k - "j77 (fkfk)MF 



*/ 



(2.36) 



In terms of these Green's functions, the stationarity equations (2.27), (2.28) are 

E Tr K,k^) + E"k%rt = M av , a^N b (2.37) 

k k 

J2 m a V k9k = ~M av , aeN f (2.38) 

k 

J2 T * G * + Y, v i v k = -Z>> ( 2 - 39 ) 

k k k 

where Tr traces over Nambu indices, the minus sign is due to the fermionic statistics of monons and we have 
written the Fourier transformation for the time-independent parameter v r as v r = ^Y i>k. exp(ikr) (here the 
sum over k does not include the factor of 1/V). The factor e , S = +0, is assumed to be inserted under 
the sums over k in these equations [10]. Now the obvious identity J2 k gf-lng^ = and Eqs. (2.32), (2.36), 
(2.38) yield 

it v (M 4v M 6v - M 3v M 5v ) = (2.40) 

or lm(M4 V M6 V ) = 0, so the phases of M^ n and Me v are opposite. Consequently, in the states with zero values 
of the parameters A v = M' ir] by the gauge transformation a vx — > a vx + Sa v with a constant Sa v we can make 
all parameters M ar) , a = 3, . . .,6 real and positive. For simplicity, we assume that they are also real and 
positive in the case of nonzero A v . We will use for these parameters the notations B v = Mb v = Ms v = M^ n 
and F v = Mf v = M<; ir] = Me v . Using Eqs. (2.38) and (2.43) below, one can easily see that F v is different from 
zero if and only if B v is different from zero [11]. 

Now the spinon Green's function (2.35) may be written in the explicit form: 

G k = [-iuJT 3 - ij, b -2t v F v cosk v + J v (A' v t 2 + A"Ti)sink v ] 

iiJT 3 — [lb — Itr/Fri COS &jj ~ J-q (^ T 2 + -4" r l) sm &!) ' 

oj 2 + (fib + 2t v F v cosk v ) 2 — \J V A V smk v \ 2 
where A' and A 1 ' are the real and imaginary parts of A v , respectively. The dispersion law of spinons is given 

by 

wjj. = \J (nb + 2t v F v cosk v ) 2 - \J V A V sin^p. (2.42) 

The monon Green's function (2.36) may be written explicitly as 

g k = (-iw - [if - 2t, 1 B rl cos k^)' = (-iw + e^.) , (2.43) 

where the quantity ej^. defined by this equation gives the monon dispersion law. 

The parameter i> r is determined by requiring the linear part of the action Su„ = J jCu„(x) dx to vanish. 
We have 

Sun = E( 6 k G k lu k + 4 G k 1& k)- (2-44) 

k 

where we have denoted G7 = G^ \ k -/ V), an d similarly for &K-- Now i>k. is to be found from the equation 
G7 v-^ = 0, which has a solution only for k = k c such that ui^ = 0. Therefore, for k ^ k c , ^ = 0, while for 
k = k c we obtain 

/ -iJ„Ar,smk c „ \ 

«k.=«k. _ p . > ( 2 - 45 ) 

\ Hb + tt n r n cosk crl J 

where a^ is a coefficient of proportionality. Evidently, if wj^. = then also wj^. = 0; the coefficients a+k 
along with the other parameters are to be obtained from the self-consistency equations (2.37) - (2.39). We 
note that these equations determine only the sum |ai- | 2 + |a_i c | 2 - Therefore, with this sum given, the phases 
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of a+k and their relative absolute values may be chosen arbitrarily, which reflects the existence of different 
possible ways of breaking the symmetry. 

Instead of performing a shift of spinon fields, one can use a different, and often more convenient way to 
describe the effect of spinon condensation [12, 13]. Namely, since we consider a large, but finite lattice, the 
expectation value of b x may be different from zero only in the presence of the infinitesimally small ordering 
external field, which points the direction of the symmetry breaking. In the absence of such field, the expectation 
value of b x will vanish, and the ordering will manifest itself in the appearance of the abnormally large (of order 
f3V) values of the boson Green's function at momentums k c = (0,k c ). Indeed, in this case the MF order 
parameters are such that uu is not exactly equal to zero, but has (for nonzero temperature) a small but finite 

value of order (j3V)~ 1 ' 2 [12]. From (2.41) we have then G kc = f3Vv-\. vi , where v-i. is given by (2.45) with 
ai. = [—/3Vui? (fif, + 2t v F v cos k crl )~\~ l ' 2 . In this way, the theory of the phase with v r ^ may be presented 
in the same form as that for v r = 0, and one only has to carefully separate out the terms with k = k c in the 
sums over k. Formally, this separation is achieved by the substitution 

G k ^G k +f3Vj2 v kA <W (2-46) 

K 

Finally, in the Fourier representation, the interaction action in (2.22) is 

& = -^^i,-X, g - (2-47) 

q a 

Therefore, the AF-monon and AF-spinon interaction vertices U ar} (k, q) have the form 

U av (k,q) = e'^m avk , (2.48) 



where for a = 7 the exponential factor in this formula must be omitted, and by m 7 \. we mean m '/ = ± 
depending on whether AF <f>^ q interacts with a spinon or a monon. 



II. 3 Gauge fixing and the AF sector 

Here we want to discuss the AF sector of the diagram technique. First we need to choose a GFT. We define 
a "left" lattice derivative d% , and its staggered counterpart d% , by the relations d^ip x = l~ 1 (ip x — ip x -fj) and 
d^ = frdfer, so that dfa x = l~ l {i) x + ^-^)- We then use the GFT 



J^, a \X) — ^.±j„. 



-g\ 



2C 
1 



(2.49) 



^x — q <J r\ \"^V I r\ y^rjx ®>r} ) T £^ D^ r ^ O^ \Q>r}x 



*VJ' 



8 

where (* is an arbitrary positive number and a v , a v are the gauge fields at the saddle point: 



a 

All t-^ 



(2.50) 

^ Oif) • 



Evidently, in the case when parameters M arj , a = 3,..., 6 are real, a v is equal to zero. The reason for 
the choice of the GFT in the form (2.49) will become clear in the next section. Although a vx and a vx are 
logarithmic functions of the original AF, we will see below that in the framework of the 1/N expansion one can 
develop a diagram technique, using the GFT in the form (2.49). It may be easily seen that the Faddeev-Popov 
determinant, corresponding to the GFT (2.49), is independent of AF and hence may be ignored. 
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The GFT (2.49) really fixes the gauge only in the phases, where for every direction r\ at least one of the 
quantities A v and B V F V is different from zero. This condition is not satisfied in the symmetric phase, which 
has vanishing A v , B v , and F v for all 77; however, in this work we will not be interested in this phase. Also, 
this condition is violated in the physically important phases with A Vo = B Vo = F Vo = for only one direction 
rjo- In this case we will consider the problem as two-dimensional, neglecting the existence of the direction rjo] 
in the plane, orthogonal to 770, our condition is satisfied, and we can use the GFT (2.49). 

We note that the GFT (2.49) does not eliminate the gauge freedom completely: there still remains the 
residual symmetry with respect to the gauge transformations, which depend only on time and are space- 
independent. However, for nonzero e our action is not invariant with respect to these transformations. There- 
fore, the corresponding divergences will produce contributions proportional to 1/e; as was already pointed 
out, these contributions must cancel each other in every order of the 1/N expansion (see Ref. [9], where such 
cancellation was demonstrated explicitly in the first order in 1/N). A different way to eliminate this residual 
gauge freedom will be discussed at the end of the next section. 

We now want to learn how to work with the GFT (2.49). For that we note that d vx contributes to (2.49) 
only in the phases with <f>i y 2ri 7^ 0, and a vx - only in the phases with <f> ari ^ 0, a = 3, . . . , 6. We can then 
use the fact that the parametrization (2.15), which is not invertible at (f) avx = 0, is invertible about the point 
(f) arjx = <j) art . The inversion is given by the expansion of (2.17) - (2.18) around <f> ari in powers of (f) avx as 



^r]x 



7p]^2go.ln(<f>a V +<f>a V x) = a v + ij> vx + ij>' vx , 

1 ^ , ( 2 - 51 ) 

-r7j22gaM < f ) ar 1 +(f>ar 1X ) = a^ + (f)^ + (f)'^, 



where 



a 



J OLT}X 

(par] 
^ar]x 
(par] 



(2.52) 



and (f)' , <f>' include all the terms of second and higher powers in AF, coming from the expansion of the 
logarithms in (2.51). 

We have now for the GFT C g (x) = C g (x) + C Jx), where 



4 0) (*) = £4 0)2 . 

Li 0) = ^\A v \ 2 d^ vx +t 2 v B v F v d^ vx , 



(2.53) 



and C'(x) includes the contribution of <f>' and <f>'„ x . We will regard C'(x) as describing the self-interaction of 
AF. Although it has an infinite number of terms, it will be clear that only a finite number of them contribute 
to every order of the 1/N expansion. 

The effective free action of AF (the sum of all terms quadratic in AF and proportional to N) is given 
now by the bare part in (2.22), by the quadratic part of GFT and by the contribution of the random-phase 
approximation (RPA) polarization bubbles: 



eSS = f[c ( °Ux) + C^(x)] dx 

■— I axax y ^ <p aTIX YY rjarj i 



2 

H-riari'a'( X ~ X ) = i^arjx^a'r] 1 x') M F 

11 



(x — x')d) n i„' x ' , \ ■ I 



With the presence of the gauge-fixing part, this action is nonsingular and unambiguously determines the bare 
propagator of AF D vav i a i(x — x') by the inversion of the corresponding quadratic form. Like the parameters 
t>r, M arj , the propagator D is determined by this inversion on the MF level and is different from the full 
dressed Green's function ((/) aV x4 l a l r] l x 1 )- Since we have used RPA bubbles as a part of the free action of AF, 
the diagrams which include the segment shown on Fig. 1(b) must be excluded. 

Obviously, the propagator D is proportional to 1/N , so the order of every diagram in the 1/N expansion 
is equal to the number of propagators D minus the number of b- and /-loops and £' -vertices in the diagram. 
Also, if we integrate out the b- and /-fields, the resulting effective action of AF will be proportional to 
N. Consequently, the order of any diagram in the 1/N expansion is equal to the number of AF-loops in 
this diagram. The MF free energy per color is of order 1, and the diagrams give to it corrections of order 
(1/N) P , p > 1. Therefore, the MF contribution is the only one that survives the limit N — > oo. As for the 
corrections, the simple analysis shows that there exists only a finite number of diagrams contributing to each 
order in 1/N , which makes it computable at least in principle. 

We now write the gauge-fixing part of Eq. (2.54) in the form 



-J d^'X^«^ n w«<( 



4°)(.) dx=-J dx dx > ^ ^IC^(* - *')*«•**• > (2-55) 



which determines the matrix IT?^ „' a i{ x )- Then in momentum representation, we have for the effective action 
of AF (2.54) 

^eff = 71 2^/ L^t ^arj-q^rjarj'a 1 (?) <t>a'r)'q , (2.56) 



2 

q a,a' 



where the polarization matrix is 



iW«< (g) = nj°2 a , <*„,,< + nj5 a ,(«) - n' vav , al (g) 

^a^'a' (q) = -ay(M a r,q M a'ri'-q)MF, 



(2.57) 



and H 9 na Ji a i(l) is given by the Fourier transformation olH 9 a l,(x). In explicit form, the nonzero elements of 
LV in the phase with v r = are 



vav'a' W/ 



' #e ^-V)]rTr (m a „ k G fc+f m aVk G fc _ f ) , a,a'etf b 



k 



where in the top (bottom) line of this formula m^ is defined as mK(m , ) and LV 77 (g) is equal to the sum of 
expressions, given by the top and bottom lines. The propagator of AF is 

V«'(?) = [n(«)C,v • ( 2 - 59 ) 

As was explained above, in the absence of the ordering field the same formulas (2.58) describe nonzero elements 
of LT in the phase with condensed spinons, where the substitution (2.46) has to be made if the explicit separation 
of condensate contribution is desired. 

III. Collective modes 

The breaking of continuous symmetries implies the existence of collective modes, which we will study in this 
section. We will use here the formulation without the ordering field, which unifies the cases of zero and nonzero 
i> r . The initial theory has the gauge invariance with respect to the transformations (2.19) - (2.21) with an 
arbitrary function \x ■ We want first to consider the phase with A v = and nonzero B v , F v . After fixing the 
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gauge, the theory is invariant with respect to the residual symmetry, described by Eqs. (2.19) - (2.21) with 
X x = Xv x v The corresponding transformation of a vx is a global shift a vx — > a vx + l~ 1 Xv> so that AF <f> avx , 
a = 3, . . . , 6 undergo the global phase transformation (f) avx — > exp(ig a Xri) < f>arix- We have therefore the global 
symmetry, described by a continuous group with d parameters \n ■ The nonzero values of B v , F v correspond 
to spontaneous breaking of this symmetry. Let us show that this symmetry breaking leads to the generation 
of d massless Goldstone modes, which are just the Cartesian components of the corresponding "photon" a vx 
(compare with Ref. [14]). 

In the phases with a broken symmetry we denote by <f> a 'ri'q the deviations of AF in momentum representation 

from their vacuum values f3V<f> a ' v '3 q o- We want to find the values 4> a i''J , which correspond to the infinitesimally 
small field a vx (where r\ is kept fixed at some particular value) with zero momentum. These values are the 
deviations, which result from the shift of a vx by the infinitesimally small value Sxrj from its vacuum value, if 
all other parameters of AF are fixed at their vacuum values. From Eq. (2.15) we have then 



^a'l'O — ^V'v9a'<f>a'ri', 



(3.1) 



where we have dropped nonessential overall factors. 

We must now prove that this mode is an eigenmode of the effective Hamiltonian of AF with zero eigenenergy. 
This means that <f>^l must satisfy the equation J2 a , H v " a "ri' a' {Q)4 '„/„/ = ^> wl th ^-ri"a"ri'a'{<l) given by Eq. 
(2.57). Since in this phase nf,,^„ /Q ,,(g) does not contribute for q = 0, we must show that 

/_^ / N-ri"a"ri'a'\Q) < i > a<'}i'0 = /_^ ^ri"a"a' ^«V'0 • (3-2) 

a' a' 

To do this, we consider the theory with the Lagrangian 

C K {x) = Cmf(x) - ^2K a i v i x M a i v i x , (3.3) 

a' 

where K a i v * x are the sources for M a i v i x . Let Zk be the partition function corresponding to this Lagrangian: 



Zk = V(b*bf*f)exp - / £ K (x)dx 



(3.4) 



We now make in the path integral in (3.4) the change of integration variables (2.19) with \x = Xr\ x n (since 
r\ is fixed, there are no summations over r\ here and below!). We have then the transformation law for the 
composition operators M a i v i x r. 



M, 



a'r/' x' 



exp {ixn [{la 1 (2^ + 8^^) - ga'^'n ] } M a i^ x i 



(3.5) 



We recall that A v = 0, so in £mf(%) there are no terms with Mi^'x'- 

The condition that Zk does not change with the change of integration variables yields then 







SIiiZk 



Xv = ° 



H-iXv 

Z K X I dx / V(b*bf*f)< y^ijla'^gc'Ma^v^^ 

+ X] K a' v 'x'M a > v > x > [g a >S v i v -g a i(2x' v +S V > V )] > exp 



(3.6) 



Ck{x) dx 



The first term in the integrand will not change if we write there M' a , , x , instead of M a i v i x i . We then differentiate 
(3.6) over K a ii v ii x n and set K avx = to get 



2_^^-'r)"a"r) l a'{fy&ri'riga' < i>a'ri' + $r]" rj N M a " v " g a " — 0, 



(3.7) 
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where we have used the self-consistensy equations and the zero value of A n in our phase. Rewriting the second 

term in (3.7) as S v » v J2 a < 9a"'tt ri J a H a i<f>a<ri" = -<V'jj J2a< ^-n''a"a'9a'<f>a'ri" we obtain the desired relation (3.2). 
We note the crucial importance for our proof of the absence of the terms with Mi^x in Cmf- 

Consider now the phase with B v = F v = and A v 7= 0. Here the residual symmetry is described 
by Eqs. (2.19) - (2.21) with Xx = £rXij#jj- A nonzero value of A v means that this residual symmetry is 
spontaneously broken; we expect that the d massless Goldstone modes generated by this symmetry breaking 
are the components of the "staggered photon" d vx . The form of the residual symmetry suggests that in this 
case the components of d vx with momentum q = 7? have zero energy. To check this assumption we first obtain 
from Eq. (2.15) that the values <f>^l x corresponding to infinitesimally small field a vx with momentum 7? are 

proportional to S ri i ri £-r9a' < f>a'ri' ■ One can show now that the corresponding amplitudes <f> ,,-. = S v i v g a i<f) a i v i 
satisfy the equation 

$>;» a - W (0,7T)^ = X) n S»a»a^S»ff. ( 3 " 8 ) 

a' a' 

which may be proven by the same method as (3.2) using the function Xx = £.rXv x v m (2-19). Since in this 
phase the GFT does not contribute at q = 7? , this means that the field d vx with momentum 7? has zero energy. 
It is again crucially important that in this phase Cmf does not have terms with M arjx , a = 3, . . . , 6. 

Apriory we cannot exclude the possibility of the existence of more general "mixed" phases, which for some 
77 (we will denote them rj u ) have A v = and nonzero B v , F v , while for other 77 (which we will denote rj s ) have 
A v 7= and B v = F v = 0. To describe these phases, it is convenient to introduce parameters i> v , v n such that 
i> Vu = D Vs = 1 and v r)s = v rlu = 0. The corresponding Goldstone modes now are a VuX and d VsX ; they have 
zero energy in the state with momentum k m , where k™ = itv^. To show that, one should repeat the same 
consideration as above, using the function Xx = exp(ik m r)x^a;^. Then one obtains 

^n;„ a , Va , ( o,k™)^) km = £<U<V> (3-9) 

a' a' 

with 4', V ,\ cm = $ri'ri < f>ct'ri'{vri'9ct' + ^rj'Qa 1 ), which proves our statement. 

Now we want to consider the phase with nonzero A v and B v , F v . In this phase the situation will be 
quite different. Indeed, the Lagrangian Cmf in this case includes the terms with M ar]x for all a = 1, . . .,6. 
One can easily see that this results in a failure of our consideration above: Eqs. (3.2), (3.8), and (3.9) will 
not be satisfied. This is because in this phase the fields d vx and a vx got screened, obtaining a gap due to 
the Anderson-Higgs mechanism. Indeed, consider first the phase with A v = and nonzero B v , F v . As we 
have just shown, in this phase a vx is a gapless gauge field. The state of the gauge field with zero energy 
corresponds to such a momentum, that the gauge transformation with this momentum does not change the 
gauge field. For the field a vx this momentum is equal to zero, and the corresponding function Xx in (2.19) 
is a constant: Xx = X- On the other hand, under the gauge transformation with this function Xx, the AF 
<f>i,2r)x acquire a factor exp(— 2igi2x)- Therefore, they will produce a gap in the spectrum of the gauge field 
a vx , if we cross a border and pass into the phase with nonzero expectation value of <j>i t 2r)x, i- e - with A v 7^ 0. 
Similary, in the phase with B v = F v = and A v 7^ the gapless gauge field is d vx . Its zero energy state is the 
one with momentum 7?, since it does not change under the gauge transfromation with a function Xx = X^r- 
Here, the fields 4> aV x, ol = 3, . . .,6 are not invariant with respect to this transformation, acquiring a factor 
exp(— 2i^g a x)- Consequently, in a phase where they have nonzero expectation values, they will give a gap to 
the state of d vx with momentum 7?. 

When the polarization matrix 11(g) has an eigenmode with zero eigenvalue, the inversed quantity - AF 
propagator D(q) - has a pole. As is well known [15], the poles of the AF propagator correspond to the bound 
states in the original system of interacting particles, in our case - spinons and monons. In the case of d vx this 
bound state exists in the particle-particle boson channel, while in the case of a vx it exists in the particle-hole 
channel of both bosons and fermions. It is important to note, however, that although we found d eigenmodes, 
only d— 1 linear combinations of them correspond to genuine Goldstone modes. Indeed, one linear combination 
of these modes corresponds to the possibility of performing the gauge transformations. This mode appears in 
all the phases with a broken gauge symmetry; in the polarization matrix II' it is completely decoupled from 
the gauge-invariant operators and, in the absence of GFT, has zero frequency for all momenta [16]. 
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To find the explicit expression for this zero mode in momentum representation we make in the path- 
integral (3.4) the change of integration variables (2.19) with the function Xx = X ex -P(iqv x v)> wnere this time 
the summation over r\ is assumed. For small x, the variations of different terms of Ck are given by 

8M a , v , x , = x le igr ' Xr > (d qn rg a , + d qv ,g a ^ M a , n , x ,, (3.10) 

where d qv and d qv are the Fourier images of the usual and staggered lattice derivatives: 



1 (l _ e W') 



i = -d 



e 



■q v ] 



2 ,-ii . q v 

— e 2 sin — - 

/ 2 q v - 


, In 
+o' I 


2i ill q v 

— e 2 cos — -. 

/ 2 





(3.11) 



Therefore, we obtain by the same method 



E 



6 (o) 

y a'i)'q 



{d qv ig a i + d qv ig a ij (f) a i v i . (3.13) 



do) 



These equations show that for every q in the absence of GFT the mode <j> a ln is an eigenmode with zero 
eigenfrequency. A different derivation of Eq. (3.12), using the properties of the b and / Green's functions, will 
be presented in Appendix. 

Equations (3.12) - (3.13) can be also used for the alternative derivation of the results, discussed above. 
Indeed, if for every r\ only <f>i y 2ri or <f>ar), a = 3, . . .,6 are different from zero, then we can obtain (3.9) from 
(3.12) by differentiating (3.12) over q v and setting q = k m . If, on the other hand, at least for one r\ both 4>i,2r] 
and 4> arl , a = 3, . . . , 6 are nonzero, then this derivation does not work, and the only mode with zero frequency 

is the zero mode ^ q . 

Now let 1Z be an arbitrary gauge-invariant operator. By adding the term K 1Z to the Lagrangian Cmf and 
considering the variation of Zk (3.4) with respect to x and K, it is easy to get 

<ttE^£VC,qW = ' ( 3 - 14 ) 

a 

thus proving that the zero mode and the gauge-invariant quantities are not coupled through the matrix II'. 

Now we want to make some comments on our choice of the GFT. In the real space, the states of the zero 
mode have the amplitudes <f> a 'ri'x', which are equal to the deviations of AF from their vacuum values <f> a 'ri', 
generated by the gauge transformations (2.19) - (2.21) with arbitrary functions Xx' ■ For nonzero e, these 
functions must be time-independent. A convenient basis set of them is given by the (5-functions x x ' — X^r'r- 

From Eqs. (2.15) and (2.19) - (2.21) we have for small x the expression for the corresponding zero-mode 

(r) 
amplitudes <j> a , v , x ,: 



■l>a>i> x > = il x4>a'n' [~g a 'd* + g al df, ) 8 T , T . (3.15) 

In the space of AF, taken at some particular moment of imaginary time r, the scalar product is defined as 

6 
A(l) A&)\ = V 



i) /_ j \J ar\x ■ T ar\x ' T ar\x ■ T ar\x J ' ^ix^lx 

r' L a=l 

1 V^ V^ jM) jS 2 ) 

a/3 r' 



(3.16) 



where in the right hand side of this equation x' = (r, r'), <j> a „J and <j>' x i are the real and imaginary parts 

of the AF <f> a ^J, jajS is a matrix with nonzero elements 712 = 721 = 735 = 753 = 746 = 764 = 777/2 = 1, and 
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we have used Eq. (2.12). In the absence of the GFT, the quadratic form of the effective free Lagrangian of AF 
in Eq. (2.54) does not depend on the components of AF along the zero mode states <f>( r >. These components 
are equal to the scalar product (<f>, <f>( r >) T . The lattice derivatives d^' L , d^' L allow integration by parts: 

]T f x V f g x = - Y, ^ v L n h , (3-17) 

r r 

where \J V stands for both d v and d v and the functions f x , g x satisfy the periodic boundary conditions (or the 



r 



ration 



lattice is assumed to be infinite). Now using the expression (3.15) for <f> a iL x i and performing the integ: 
by parts in the definition (3.16) of the scalar product we obtain the equality 

Therefore, the quadratic part (2.53) of the GFT gives the potential energy to the component of AF, which 
originally (i.e., without gauge fixing) did not have potential energy because of the gauge invariance. We note 
that one can also regard C g (x), Eq. (2.53), as the GFT, with C (x) being the logarithm of the corresponding 
Faddeev-Popov determinant. In the limit C, — > this leads us to the formulation used in [16, 17]. Keeping C, 
finite, however, allows us to obtain explicit expressions for all contributions, corresponding to the gauge-fixing 
procedure, which will be important for us in the forthcoming sections. We also note here, that rather then refer 
to the finite e in order to prohibit the time-dependent gauge transformations, one can eliminate the residual 
gauge freedom with respect to such transformations by adding to L x (2.49) the term (2/ 2 ) _1 9oaoj;. Then for 
the corresponding new L x , Eq. (3.18) will remain correct, if in its right hand side [cj>, <f>( r >) is changed to 
((f), (f)^" 1 ) , where 

H^a'r,' [jja'dfi + 9 a' df, J + S a < 7 



*(*) 

->„/„/„/ 



W(r'-r) (3.19) 

is a zero-mode amplitude, generated by the time-dependent (5-function x x i = X^r'rfH - ?"' — r ) an d 






is a scalar product of AF, taken on the whole interval (0, /?) of the variation of imaginary time. 

IV. Equations of motion 

The variation of the Lagrangian (2.13) with respect to AF yields the useful equations of motion 

p 

The use of these equations, however, requires special care and in a simple form (4.1) is possible only under 
special conditions. Indeed, our diagram technique is produced by the Lagrangian (2.22), rather than (2.13). 
To elucidate the possibility of using Eq. (4.1) in our diagram technique, we start with the obvious equality 

V(b*bf*f </>)—?— e-J £{x)dx = 0. (4.2) 

With the Lagrangian (2.22), Eq. (4.2) leads to the identities 

«„) = (EO- + #^)' ( 4 -3) 

r, T (XT) X 

16 



which are just the manifestation of the equations of motion in our model. In the last term in (4.3) S g = 
J C g (x)dx is the gauge-fixing part of the action. Here we want to show how our diagram technique meets 
(4.3). We will employ the generalization of the method, used in [9]. We start with the observation that the 
operator M' in the left hand side of (4.3) can produce three kinds of (sub)diagrams, presented in Fig. 2(a-c). 
On the diagram of Fig. 2(a) the b-f lines, leaving M' , go to different vertices. On the diagram of Fig. 2(b) 
these lines, after leaving M' , meet again in the vertex of interaction with AF, while the AF then propagates 
till the next vertex of interaction with the b-f line. Since this last vertex is produced by the action of Ci n t, 
and not by an application of the external field, the b-f lines, leaving this vertex, must then go to different 
vertices. Physically, the first diagram describes the bare action of M' x , while the second - its renormalization 
due to the screening. Finally, on the diagram of Fig. 2(c) AF propagates till the vertex produced by C'{x). 



The analytic expression corresponding to this vertex is — dS' /d<j> a "ri"x" , where S' 
is the label of the AF line entering this vertex. 



M' 

arjx 



M' 

arjx 



)dx and a"r\"x" 



M' 

arjx 

•— — 





(a) 



(b) 



(c) 



^-^rjarj'a 1 \% X 



L± T]aT]'a'\' 1 ' "° / 




5S> 



6<t> 



arjx 




(d) 



(e) 



(f) 



FIG. 2. Diagrammatic representation of the equations of motion, Eq. (4.3) 



Now introducing the "bare" polarization matrix 

t(0) , 



rt v b L a ,(x - x') = n[ l,s vv ,s(x - x>) + Ki%,(x - x>) 



and using the relation 



^- V ari'a'{ x — x) — II „, a ,(a; — X ) — (D ) , a , (x — x ) 



(4.4) 
(4.5) 



17 



between the polarization matrices and the AF propagator, we can see that the sum of the diagrams Fig. 
2(a-c) is equal to the sum of the diagrams Fig. 2(d-f). Since this latter sum represents the action of the right 
hand side of Eq. (4.3), we conclude that the calculations in our diagram technique automatically satisfy the 
equations of motion (4.3). 

By the same logic one can prove the more general identity 

where T is an arbitrary functional of AF and the fields b* , b, f* , /. In the diagrammatic proof of Eq. (4.6) 
the diagram of Fig. 3, which exists when T depends on AF, should be added to the diagrams of Fig. 2(a-c). 
It is important that since all the diagrams of Fig. 2 are of the same order in 1/N, the identities (4.3), (4.6) 
are satisfied in each order of the 1/N expansion separately. In the special case of a = 7 and T independent 
of (f>7 X , (4.6) gives N~ l (TM*i x ) = 0(e); the factor N~ l here cancels the summation over colors in M-? x . We 
conclude that in the limit e — > our diagram technique satisfies the constraint (2.10) as an operator identity 
in each order of the 1/N expansion separately. 




M' 

arjx 



FIG. 3. The diagram for (T M' ar]x ), which should be added to the diagrams of Fig. 2(a-c) in the case when T depends 
on AF. The shaded circle represents the diagrams for Sy 7 /8(f> a i v i x i , where aiji is the label of the AF line entering the 
circle. 



Now we want to show that in many important cases the term with 8S g /6(f) avx in (4.6) gives vanishing 
contribution. We first show that L x , Eq. (2.49), is decoupled from the gauge-invariant quantities: 

(L X TZ) = (4.7) 

for arbitrary gauge-invariant functional 1Z. To prove this identity, we add to the action the term K1Z, where 
K is a source of TZ, and make in the path-integral representation for the partition function the change of 
integration variables, corresponding to the gauge transformation with the infinitesimally small parameter Xx- 
In the state with a broken symmetry this change is given by the formulas 

bax — > b ax + i\x(var + box), 

<t>a^x > <t>a n x+il(<t>a n + <t>a n x) (i/adf + QadfjXx , «=1,...,6 

and the old formulas (2.19) - (2.20) for the transformation of f x and 4>7 X . The variation of the gauge-fixing 
part of the action under this transformation is 



N f 
SS g = — dxL x A\ x , 



(4.9) 



Since the rest of the action is gauge-invariant, the condition that partition function does not change with a 
change of integration variables yields 

= SZ = - I V(b*bf* f<f>)SS g exp(-5 + KTZ), (4.10) 
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where SS g is given by (4.9). Now differentiating this equality over A\ x and K, and setting K = 0, we obtain 
(4.7). Along with (L X 1Z), the derivatives of (L X 1Z) with respect to x also vanish. 

Now returning to Eq. (4.6), we have there after integration by parts and using Eq. (2.51) 



SS g N J^\A v \ 2 g a d^ + 4t$B ri F ri g a d? 



Ha^x 16i/C 



L x , (4.11) 



ar]x ±\jii^ Yocr] ~ Yocrjx 



where the term (<t> ari + <f> ari x) 1 is understood as a power series. Under the gauge transformation (4.8), the 
term <f> av + <f> aV x in (4-11) obtains the factor exp[il(g a d^ + g a d^)xx]- Therefore, if the functional T in (4.6) 
obtains the same factor, then as was just shown, the term T8S g /8<f) artx there may be neglected. For example, 
it will be so in the important case when T = Mp nx with /? such that h na p ^ 0. In such cases, instead of Eq. 
(4.6) we can use a simpler relation 

(TM' avx ) = <V;c n iV^) - {^-)- ( 4 - 12 ) 

This simplification, allowing us to neglect the contribution of 8S g /6(f) avx in Eq. (4.6), is a usefull property of 
our GFT. Below, when we use Eq. (4.12), it is always assumed that the conditions for neglecting the term 
with 8Sg/6(f) avx in (4.6) are satisfied, so we will omit the corresponding discussion. When 8T /S(f) avx = 0, Eq. 
(4.12) allows us to use the equations of motion (4.1), following from the Lagrangian in (2.13), as operator 
identities. The meaning of this rule is that the total (screened) action of M' comes from the creation of the 
intermediate state of the b-f plasmon. The plasmon is represented by the set of AF, and the amplitude of 
transformation of M 1 to <f>p vx is equal to ITlJg. 

V. Interaction with electromagnetic field 

We now want to introduce to the model the interaction with the external classical electromagnetic field. The 
interaction with the scalar potential (p(r) is described by adding to the Hamiltonian (2.1) the term 

e J2 <p(r)Ar^r = \ J2 ?(*) (4&r + flfr) , (5.1) 

r r 

where e is the electron's charge and we have used the constraint to make the transformation. The generalization 
1 — > N does not affect this form of interaction. The interaction with the vector potential A(r) is obtained in 
a standard way [5] by modifying the hopping term H t in (2.5) as 

H * -> "I E [e lelAr ' ir) bl + ^rfl + ,fr + e-^-VbU^flfr+t] , (5.2) 

r 

where we made the 1 — > N generalization. The electric current is defined by the usual relation j„(r) = 
— SH/SA V (r) . Up to terms of first order in A, we have j e (r) = ej(r) +j d (r), where j(r) is the electron current, 
while ej(r) = j p (r) and j d (r) are respectively the paramagnetic and diamagnetic parts of the electric current. 
We have 

hi r ) = % -fi- (blr+fjbarfl+fjfr ~ bl r b ar +^flfr+^) , (5.3) 

J» = - ! X L A('') (i + ,-^4/r + t^r + jj/r + ,-) • (5.4) 

In terms of the operators M avr , the electron current is 

j v ( r ) = ■^'52 M av rj a h va pMp v r , (5.5) 

a/3 
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where the j a were determined in Eq. (2.16). Due to the presence of the EMF the action obtains an additional 
contribution 



dx < p(x)ip(x) 



i$(*) + ^(z) 



A,{x)\, (5.6) 



where we have neglected the terms of third and higher powers in elA v , and have denoted p(x) 
It is convenient to make in (5.5) the substitution (4.1) and obtain 



[J$(x) + Jf(x)] 



(5.7) 



2 L " v 
where we have also used (2.24) and have defined 

J^ f (x) = il^2j a M avx (<j> av +<f> avx ). (5.8) 

The currents J^ (x) defined in (5.8) are invariant with respect to the gauge transformation (2.19) - (2.21); 
one can identify them with the currents of spinons and monons. We will soon see that this identification is 
not useful, but first we will consider these currents in more detail. We note that these currents may be also 
obtained in a different way. Namely, we can modify the Lagrangian in (2.11) in the same way as in (5.2), 
and then decouple it using AF 4> avx - The factors exp(±ielA vx ) should then be decoupled to the products of 
exTp(±ieKi}lA vx ) and exp(±ieKflA vx ), where Kb (k/) is a part of the electron charge, belonging to the spinon 
(monon). All results must be independent of the choice of Kbj, if they satisfy the condition Kb + Kf = 1. The 
term Y^ a §anxM ar]x in (2.13) will then obtain the form Y^ a ex -P( ie lJaK a A vx )(f) avx M avx , where n a = K b (nf) for 
a G Mb{Mj). Now differentiating this Lagrangian with respect to — eA vx for A = 0, and performing the shift 
of AF, we obtain for the electron current 

j v (x) = K b J%(x) + Kfjf(x). (5.9) 

Equation (5.7) corresponds to the symmetric choice Kb = k/ = 1/2. 

For the current (5.9) to be independent of the choice of Kbj, we must have 

J b v (x) = jf(x) (5.10) 

as an operator identity. It is easy to see that this equality is a consequence of the absence of nontrivial bare 
kinetic energy of AF in our theory. Indeed, we can for example consider the variation of the Lagrangian (2.13) 
with respect to the transformation a vx — > a vx +Sa vx and use the equation of motion SC/Sa vx = 0. Under this 
transformation, the variations of AF in (2.13) are 8<j) artx = ilg a (/) avx Sa vx . Performing now the differentiation 
of the Lagrangian (2.13) over a vx and a shift of AF, we obtain 

= SC/Sa vx = -il^2g a (<j> av + <l> avx )M avx . (5-11) 

a 

The last expression here is equal to Jl(x) — </„(«), so we immediately obtain (5.10). This also implies j v (x) = 
J^'f (x), in agreement with [18]. If AF had a nontrivial bare kinetic energy, then there would have been an 
additional term on the right hand side of (5.11), which would violate the equality (5.10). On the other hand, we 
can show that (5.10) is a consequence of the equations of motion, considered in the previous section. Indeed, 
using (4.12) on the right hand side of (5.11), we get 

Jf(x) - J*(x) = -il^2((/) av + (/> aV x)g a n- v a l 3((l>/3 V +(/>/3 V x) + il^2g a , (5.12) 

a/3 a 

which vanishes due to (2.16) and the antisymmetry of g a Tv*„. We see that although (5.10) looks like an 
equation, which puts nontrivial restrictions on the observable quantities in our theory, it is in fact an identity, 
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which is satisfied automatically due to the equations of motion in every order of the 1/N expansion and does 
not provide any restrictions. In other words, J^ (x) are just the different forms of the electron current, and 
we will not associate them with the currents of spinons and monons. 

A convenient representation of the electron current may be obtained by the separation of different terms 
in (5.7). We substitute there M ar]x = NM arj + M' and <f> ari = J^o h r]a pMp r] and use (2.40) to get 

h(x) = ^[j b v (x)+jf(x)+j*(x)+j'(x)], (5.13) 

j b v f(x) = ilt v M Lbv J2j a M' avx , (5.14) 



■"7 



a b,S 



j*(x) = ilNY^J<*M ari <f> arix , (5.15) 

a 

j 6 v (x) = il^2j a M' arix <l> arix . (5.16) 



The formulas (5.14) may be written in the familiar form 



4W = ■^(tfK*-K*-K*-tfb (rx ) 



2m* 



HW = -r-ritffZ-U-n-tfu) 



(5.17) 



1 / v 1 
2m 



where the substitution b ax — > i> CTr + b ax should be made in the equation for j b (x) if v r ^ 0, and where we have 
denoted 

<' = Sm^- ^ 

Therefore, the hopping matrix element of the monon (spinon) determines the effective mass of the spinon 
(monon); zero hopping corresponds to an infinite effective mass. 

The factor 1/2 in (5.13) corresponds to the symmetric choice m, = Kf = 1/2. If we had made a different 
choice, we wouldn't have this factor in (5.13), but instead every term in the sums over a in (5.14) - (5.16) 
would have obtained a factor n a . As we discussed above, these changes would not affect any physical results. 

By the equations of motion (4.3), the current jt(x) is different from j b (x) + jf (x) only by the contribution 
of the GFT. As is clear from Eqs. (4.11) and (4.7) this contribution is always of higher order in 1/N and may 
be neglected on the MF level, N — > oo. Also of higher order is the contribution of jl{x). We conclude that on 
the MF level the expression for the total electron current is 



J 



■MF 

1 



(x)=j b Jx)+jUx), (5.19) 



which corresponds to the independent interaction of EMF with the currents of spinons and monons (5.17). 
One can easily convince himself that representation (5.19) remains correct for every choice of Kbj- 

Now let us consider the interaction with EMF, which is produced by the density p(x) and the currents 
jt'f'^{x), Eqs. (5.14) - (5.15). We define the corresponding effective vertex of interaction of the b-f line 
with EMF, represented by the diagram of Fig. 4(a), in such a way that the disconnection of EMF from this 
vertex would leave us with a legitimate diagram. Two obvious terms of the effective vertex correspond to the 
diagrams of Fig. 4(b) (representing the "direct" contribution of p(x) and j b '^ (x)) and Fig. 4(c) (representing 
the contribution of jfi(x)). However, apart from these terms, the effective vertex must also include the 
"screening" term, represented by the diagram of Fig. 4(d), which describes the screened contribution of p(x) 
and j '* (x) . Indeed, the disconnection of EMF from the diagram of Fig. 4(d) would leave us with the forbidden 
subdiagram of Fig. 1(a). We note, that although the diagram of Fig. 1(b) is forbidden, the diagrams of Fig. 
5(a,b) are not, and should be taken into account. The diagrams of Fig. 5(c,d) will be discussed later. 

Because of the presence of the GFT, we must also introduce the effective vertex of interaction of EMF 
with AF. This vertex is presented in the Fig. 6. We will denote this vertex as V/j ■ (<?), where its dependence 
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(a) 



(b) 



(c) 



(d) 



FIG. 4. The effective vertex of interaction of the b-f line with the EMF, corresponding to the currents j v (x), Eqs 
(5.14 - 5.15). The dashed tail represents the EMF. 




(a) 






(c) 



(d) 



FIG. 5. The (sub)diagrams which become forbidden ones after the disconnection of the EMF line. 
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on momenta other than the momentum q of the EMF is made implicit and the number n of AF, leaving this 
vertex on the left hand side diagram of Fig. 6, must be greater than or equal to 2. The n + 1 — particle 
interaction vertex on the right hand side diagrams of Fig. 6 is produced by the interaction part of the GFT, 



C' g (x) 






n n n 

FIG. 6. The effective vertex of interaction of the AF with the EMF; n is the number of the AF lines leaving the vertex. 



The analytic expression corresponding to the segment of Fig. 4(a) is 



e^, [p(-q) V bJ (q, k) - A v (-q) V^ (q, k)] F, 



fc + f ' 



(5.20) 



where F b '^ are the corresponding Green's functions: F b = G p and F^ = g p . The inte- 
raction vertices here are sums of the contributions of the diagrams of Fig. 4(b - d), 
V^'f(q,k) = ^2i Vfi (q , k) , where the index i = b,c,d corresponds to the diagrams of Fig. 4. By Eqs. 
(5.6), (5.13) - (5.15), and (2.31) we have 



v b - m ( q ,k) 



1 

2 

V bJ{b) (q,k) = t -t ri M fibri e i ^J2^ m a^k = ^ ei ^J2 ja( f > ^ m ^k' 



a b,S 



a b,S 



V b /^(q,k) = (l-S^-'-Nj^JaMa.W^iq^) 



(5.21) 



V, 



bj(d), 



k ) = N J2J2 TT [K {b) (9,k , )G kl+ ^e-^m arikl G kl _^W^(q,k) 
k' 



OLb k' 



at k' 



Here, like in Eq. (2.58), in the top (bottom) line of the last formula m^ is defined as tjiK (m , ), and 

W% (q, k) = Y, D„ ar ,,p (q) Up v , (k, q) (5.22) 

corresponds to the propagation of AF and its interaction with the b-f line on the diagrams Fig. 4(c,d). 
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Consider finally the interaction with EMF, produced by the current jl{x), Eq. (5.16). In j s , "s" stands 
for "symmetric". Unlike other terms in (5.13), this current exists in both symmetric and asymmetric phases, 
so that in the symmetric phase with M ar) = only this term gives a nonvanishing contribution to j v (x). 
The current jl{x) creates the vertex with a vector potential, inserted at the 4>-bf interaction vertex. The 
corresponding segment is presented on Fig. 7. It corresponds to the analytic expression 

-eAr, (-q) F b k [^ Y^ Va i ( k >P> ( i) ^W"' (~P) F k+s±z> ( 5 - 23 ) 

where the functions F b '^ have the same meaning as in Eq. (5.20), and 

// • In -\~Prt 

V av (k,p,q) = -j a m arik e 1 ^^. (5.24) 

We note that just like with the diagrams of Fig. 5(a,b), the diagrams of Fig. 5(c,d) are legitimate, although 
the diagram of Fig. 1(b) is forbidden. 



b,f ! b,f 

— > 

K ~r 2 ? K 2 



IP 

FIG. 7. The vertex of interaction with the EMF, corresponding to the current j^(x), Eq. (5.16). 

VI. Ward identities 

In this section we will show that the amplitudes of interaction with EMF in our theory satisfy the Ward 
identities 

iuiT (q, k) + d_ qn T n (q, k) = t 3 G^' +3 _ - G^K r 3 . (6.1) 

Here ui is the zero component of q = (w,q), and we have introduced the amplitudes corresponding to the 
electron's interaction with EMF 



f3V 

pv 



T (q,k) = — <n,*«_ f ** t + l f ) J (6.2) 

T v (q,k) = ^<j,(g) *£**£'*>, (6.3) 

r,ec 



where £ and ( are color indices, ^| is a Fourier component of the electron field 



*■($:)■ 



(e) 

and G k is the electron Green's function: 



?( e ) = J_/,T,fC,T,fCt\ 



Gr = -5^(*«** 4 '>. (6.5) 
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The summation over repeating color indices is assumed. 

In the rest of this section, we present the path-integral and the diagrammatic derivations of the Ward 
identities (6.1). It will be clear from these derivations that the Ward identities are satisfied in each order of 
the 1/N expansion separately. Only the Ward identities themselves, and not these derivations, will be used in 
the forthcoming sections. 



VI. 1 Path-integral derivation 

We start by introducing the source terms and considering the theory determined by the Lagrangian 

C K (x) = C(x) - A'«t xpCC _ *|Ct #JC ; ( 6 . 6 ) 

where C(x) is the Lagrangian in (2.11), \?|f is an electron field (6.4), and K^S is its source: 



K t< 



K 
K 






(6.7) 



The summation over repeating color indices is assumed; in the following, we will make these indices implicit. 
Let Zk be the partition function, corresponding to the Lagrangian (6.6): 



Z. 



K 



V(c)e- S " 



(6.8) 



where Sk = J £k (#) dx and c x = b x ,6* , f x , /* is the general notation for the fields, participating in (6.8). 
We now make in the path integral (6.8) the change of integration variables b x — > e la * T3 b x , b* x — > b x e~ la * T3 , 
fx -» e 8 "* f x , fx -^ e ~ la * fx ■ We nave f° r small a x the variations Sb x = ia x r 3 b x , Sb x = ia x r 3 b x + ia x t 3 b x , 
etc. The condition that Zk does not change with the change of integration variables yields 



Vic) — e K 

1 ' S(iax) 



0. 



(6.9) 



We now take into account that the Lagrangian C(x) depends on the products of fields at the neighbouring 
points, so it has the structure 

C{x) = C x (c x , C x , {Cx+fj , V = 1, . . . , d}) . 



(6.10) 



We have then: 



6S, 



K 



dx\J2 



dC x Sc x 



8c 



dC x 8c x 



dx i 



E 



dc x 8(ia x ) 
dC x Sc x 



Oi^x — rj 

dc x J 8{ia x ) ^ dc x S(ia x ) 
- 2ia x (k\ t 3 y x - *| T3 K x 
dC x 8c x+fj dC x 8c x 



E 



dc x 8{ia x ) dc x+fj 8(ia x+f j) dc x 8{ia x ) 
dC x -fj 8c x dC x 8c x+f j \ d f dC x 8c x 



(6.11) 



dc x 8{ia x ) dc x+fj 8(ia x+f j) I dr \dc x 8{ia x ) 



-t 



2 [K x t 3 V x - ^ x t 3 K x 



t. 



Now (6.9) gives 



Q w 

{Yx + id L n ff(x) - °-^r + *t T3 K x - K\ r 3 V x ) K = . 



(6.12) 
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Here (-)k means the average in the theory with the Lagrangian Ck {%), Y x = 2 sda * — T> wnere 



5C X _ \-^ 



^(ittconst 



dC x 8c x dC x 8c x+fj dC x Sc x 



dc x 8(ia x ) dc x+fj 8(ia x+f j) dc x 8(ia x ) 



(6.13) 



is the derivative of the Lagrangian with respect to the constant, ^-independent phase, and the current j^f (x 
and density n% are given by 

il ir-^ dC x 8c x+fj 



ifM 



w 



a y^ ul, x oc x+ fj (6 14) 

2 ^ dc x+i j 8{ia x+fj ) 

1 y^ dC x 5c x 

2 ^ 3c x S(ia x ) ■ l °- ° j 



It is straightforward to check that for the Lagrangian C from Eq. (2.11), jj/ (x) and n x defined by Eqs. (6.14) 
- (6.15) coincide correspondingly with j v (%) and n x . 

Since C is invariant with respect to the phase transformation with a constant phase, the first term in Eq. 
(6.12) vanishes, and Fourier transformation of this equation yields 

(iujn q + d_ qv j v (q) + ^ (*L§ r 3 A 'fe + § ~ A 'I-a r 3*fe + §))K = 0. (6.16) 

k 

Now differentiating this equation with respect to -?Q_g and KjSg. an d setting K x ' = K^S = we obtain the 
Ward identities (6.1). 

VI. 2 Diagrammatic derivation 

Here we want to show how our diagram technique satisfies the Ward identities. We first note that the form 
(6.1) of the Ward identities is not in a straightforward correspondence with our diagram technique. Indeed, 
as we have just seen, this form corresponds to the Lagrangian (2.11), rather then to the Lagrangian (2.22), 
which generates our diagram technique. With the Lagrangian (2.22) we obtain a different expression for the 
current (6.14) and for the Y-term (6.13). The current (6.14) with the Lagrangian (2.22) is 

U / (6.17) 

h( x ) = ^/2^ a [^M^ + <?W ( NM a V + M' arix )] . 

a 

The new contribution to the current, —j v (%), may be also obtained via differentiation of the Lagrangian with 
respect to — eA v (x), if along with the substitution (5.2) we also make in (2.22) the change of integration 
variables (f) avx — > exTp[—ielg a A v (x)/2](f) avx . By the equations of motion (4.12) this new current vanishes; 
however, in our diagrammatic proof of the Ward identities it will be easier for us to keep it in the formulas. 
The action of j v (%) is described by the same diagrams Fig. 4 - 7 as the action of j v (x). The corresponding 

analytic expressions are obtained from Eqs. (5.20) - (5.24) (the second equality should be used for V v (q, k) ) 
by the simple substitution j a — > g a ; we will mark the corresponding vertex functions with a tilde. 
The Y-term (6.13) corresponding to the Lagrangian (2.22) is 

Y ( x ) = -IJ2M x )- (6-18) 

v 

Performing the Fourier transformation of Eq. (6.12) and differentiating it with respect to the sources, we 
obtain that the amplitudes in our diagram technique must satisfy the equation 

T S G^ +i - Gjg f r 3 = iioT (q, k) + d. qv [r„ (q, k) - f „ (q, k)] - | ]T f „ (q, k), (6.19) 

v 
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where 



r,(?^) = ^(j,(?)*Lt[ ( + t). (6.20) 



The Ward identities (6.1) follow from (6.19) immediately, since T v (q, k) vanishes by the equations of motion. 
Using the equality d- qv + 2(i7) _1 = —d- qv , we now present (6.19) in a form 



(<0 n(e) 



iuT (q, k) + d_ qn T^ (q, k) + d_ qn T^ (q, k) = t 3 G\^ - G\\ r 3 , (6.21 

which is convenient for the diagrammatic proof. 

In this proof we will use the following important relation, derived in Appendix: 



iu> V C (q, k) + d_ qv V^ (q, k) + d_ qv V^ c (q, k) = < 



f \{G-^r 3 -r 3 G-l i ) ) c = b 
2 (^-§"^+§ j ' c = f 



(n) 



(6.22) 



It is evident from the derivation of (6.22) in Appendix that the vertex V/2 (q), Fig. 6, satisfies the equation 

iw Vt ] (q) + d. qv V™ ( q ) + d_ qv V^ (q) = 0. (6.23) 

Another important relation, which connects the <f>-bf interaction vertex (2.48) with the vertex (5.24), is 
d_ qv V av (k,p,q) + d- qv V av (k,p,q) 



2 [ T a U av [k + -,pj -U av [k - -,pj r«J 



(6.24) 



where f a are defined after Eq. (Al). To prove this relation, we substitute in the right hand side of (6.24) the 
definition (2.48) of U arl (k ± | ,p) and use (Al) to obtain the left hand side. 

The factors 1/2 on the right hand sides of Eqs. (6.22) and (6.24) correspond to the symmetric choice 
Kb = Kf = 1/2 in (5.7); we recall, however, that as we proved in a general form, the different choices of k^j 
will give the same results. 

We now go on directly to the derivation of the Ward identities (6.21). We will derive them in the usual 
way [19], using the fact that (with one exeption, which we will discuss later) each diagram of the electron's 
interaction with EMF may be obtained by the insertion of the EMF line into the diagram of the electron's 
Green's function. According to the previous section, the line of the scalar potential may be inserted only 
at the propagators of AF, spinons and monons, and at the vertex of self-interaction of AF, while the line 
of vector potential may be also inserted at the vertex of interaction of AF with spinons and monons. By 
the interaction vertex of EMF with spinons and monons we will understand the renormalized vertex of Fig. 
4(a), described by Vjf'f (q, k). Then the disconnection of EMF line from any diagram for the amplitude of the 
electron's interaction with EMF results in the possible diagram for the electron's Green's function. On the 
other hand, if we understood by the interaction vertex of EMF with spinons and monons only the sum of the 
vertices Fig. 4(b,c), then the disconnection of EMF line from this vertex could result in the forbidden diagram 
for the Green's function, as is clear from the segment of Fig. 4(d). 

Let us now consider a diagram for the electron Green's function (Fig. 8). In order to get the corresponding 
diagrams for the electron's interaction with EMF, the EMF line may be inserted at each propagator or vertex 
of this diagram. Let us consider first insertions at the &-line on the top of the diagram. For the evaluation 
of the contribution of each insertion to the left hand side of Eq. (6.21) we use, depending on the position of 
insertion, either equation (6.24) or the first equation (6.22). Then upon summation over all possible positions 
where the EMF line can be inserted, we see that the terms in the sum cancel pairwise; only two terms coming 
from the insertions of EMF line to the extreme left and the extreme right positions on the &-line escape the 
cancellation and give precisely a half of the right hand side of Eq. (6.21). Another half comes in the same way 
from the insertions at the /-line on the bottom of the diagram. Then, the insertions at the internal &/-loops 
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lead simply to the shift of the integration variable in the corresponding momentum integrals and give, therefore, 
zero contribution. By (6.23), the insertions at the AF propagators and at the vertices of self-interaction of AF 
also give zero contribution. Finally we note that there exists an exception from the rule of the correspondence 
between the diagrams for the electron Green's function and the electron's interaction with EMF. Namely, the 
latter diagrams may include segments of Fig. 5, and the disconnection of EMF line from these segments leads 
to the forbidden segment of Fig. 1(b). However, this exception does not spoil our proof, because the diagrams 
of Fig. 5, being the insertions at the internal loop, give zero contribution to the left hand side of Eq. (6.21), 
as was just explained. 

b 




FIG. 8. A typical diagram for the electron Green's function. Insertion of EMF line at each propagator or vertex of this 
diagram produces the possible diagram for the electron's interaction with EMF. 



VII. Currents correlation singularities 
and superconducting state 

The reason for deriving the Ward identities is that they provide us with a simple way to demonstrate gauge 
invariance of the electric current and to detect the superconducting state. In the first order in EMF, the 
Fourier component of the electric current is 

(f, (?)> = [k% (?) - Pv <W ] \> (?) - k$ (?) rb) > (7-1) 

where the electromagnetic kernels K^/ (q) are 



pv 

pv 



K V0 (?) = -oTF (JV (<l) n ~q) > 



(7.2) 



and 



Pv = N yar+fjb<rrf r+ fjfr + b ar b (jr+ fjf r f r+ fjj . (7.3) 



In terms of electron fields ^ k , the Fourier component of electron current (5.3) is 



M?) = ^^E*f-W c + § sm ^- ( 7 - 4 ) 

k 



k 
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k 

Therefore, the kernels may be expressed through amplitudes (6.2) and (6.3) as 

N 
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(7.5) 
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From the Ward identities (6.1) we have then 
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(7.7) 



where we have made a shift in the summation variable and have used (7.6). Now from (7.1) we see that our 



current is gauge-invariant: the pure gauge external field (<p(q), A n (q)) 



/, u, qv 



produces no current. 



In an external magnetic field, described in a gauge <p = by a vector potential A v (r) = A v (q) e 8( l r , the 
Fourier component of a current is given by the first two terms in (7.1), where q = (0, q). By (7.7), the kernel 
there satisfies the equation 

V A '^<(q) =Pv d qV 



(7.8) 



Now differentiating this equation over q v i yields 
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(7.9) 
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Therefore, if K n ', (q) is a regular function of q at q = 0, then (7.9) gives K n ', (0) = p n & m i , and the diamagnetic 
and paramagnetic contributions to (7.1) cancel each other near q = 0; the system is in the normal state. If, 
on the other hand, Kit, (q) is singular, so that the second term in (7.9) has a nonzero limit at q = 0, then this 
cancellation does not happen, and in the external magnetic field the current with momentum q « does not 
vanish; the system is in the superconducting state. 

The required singularity in the current-current correlation function (in the BCS theory this singularity 
has a simple form, q v qri'/q 2 ) is known to exist in the superfluid phase of the system of interacting bosons. 
Consequently, if these bosons are charged, then this system exhibits superconductivity [20]. For the case of the 
usual 3D superfluidity with nonzero expectation value of the boson field, this singularity was considered in [21]; 
the singularity in the 2D superfluid state of the Berezinskii-Kosterlitz-Thouless (BKT) type was discussed in 
[22]. 

The spinon subsystem in the t-J model may in principle exhibit both types of superfluidity. The study of 
the BKT superfluidity requires application of methods essentially different from those used in this paper [23], 
and so here we only consider in detail the case of 3D superfluidity which corresponds to v r ^ 0. Let also B v 
and F v have nonzero values, so that the vertices V v '(q,k), Eq. (5.21), do not vanish. Using Eq. (5.19) 
we obtain that on the MF level the current-current correlation function K^ e ,(q) is given by the two diagrams, 
presented in Fig. 9, both of order N (or of order 1 per color). We need to evaluate these diagramsfor q = (0, q). 
The polarization matrix II' (q) (which is present in these diagrams explicitly, as well as implicitly, through the 





FIG. 9. MF diagrams for the electromagnetic kernel K ,(q) 
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propagator of AF) includes the condensate contribution, described by the diagrams of Fig. 10. We will denote 
this contribution n' u (q). The relevant quantity is the current-current matrix element of n' u (q): 



n 



j»j h , w 



a b a' h 
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\/ 3 
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k c -q f 
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yj sin [k cn , 
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k c +q f 



k cv + y I sin 



A , + ^i 



T3t; k c 



(k c 



(7.10) 



where we have used Eq. (5.14). The Green's functions G^ , „ in Eq. (7.10) are given by Eq. (2.41), where 

we have to set k = (0, k c ± q). The denominator of Eq. (2.41) is then equal to ui'L , which is proportional 

to q 2 for small q. Now the straightforward evaluation of Eq. (7.10) shows that for small q the current-current 
matrix element of n' u (q) is proportional to q v q v > / 'q 2 '. Therefore, we can expect that the diagrams with the AF 
propagator in the intermediate state give the required singular contribution to the current-current correlation 
function. 

Ck c +q 



<^J> + 



A At 
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k c -q 



FIG. 10. Diagrammatic representation for the condensate contribution to the polarization matrix, II "(q) 



However, this does not happen for the diagrams of Fig. 9. Indeed, on the MF level our problem is essentially 
similar to the one, considered by Ioffe and Larkin [24] and by Lee and Nagaosa [25]. As was first shown by 
the functional methods in Ref. [24], when the fermion subsystem is in the normal state, the singularities in 
the diagrams of Fig. 9 must cancel out, so that the whole electron system is also in the normal state. It is 
instructive to obtain this cancellation in the framework of our approach. 

We first introduce the susceptibility matrix Knarj'a'il), corresponding to the diagrams of Fig. 9, by 



K(q)=U'(q) + U'(q)D(q)U'(q) 



(7.11) 



This matrix describes all kinetic properties of the system on the MF level. For example, the MF electromagnetic 

:urrent matrix element of K(q): 



kernel Kit, (q) is a current-current matrix element of K(q 



t(«0 



where we have used Eqs. (5.14) and (5.19). Using (2.57), (2.59), and the Fourier representation LT^ „i a i(q) = 

n, ) aL' ( W + ^-lal'a'il) °f tne bare polarization matrix 11^,^,(2; - x'), Eq. (4.4), we can present (7.11) as a 
sum over all possible chains, built up from the RPA polarization bubbles: 



v- i^-i 
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Then we separate the boson and fermion contributions to the polarization matrix II' (g) 

n'(q) = n' b (q) + n'i( q ), 



(7.13) 



(7.14) 
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where II' 6 (il'-'j is described by the boson (fermion) RPA polarization bubble. The nonzero matrix elements 
of II' 6 '^(g) are given, respectively, by the first and second lines of Eq. (2.58). Now for q = (0, q) with small q 
we can present (7.13) in a usefull form: 



A'(q) 



n'U 



l + u'U 
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n' 6 (q) 



CO 
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(7.15) 



where 



CO 

n' 6 (q) = n' 6 (q)^[n( fc )- 1 (q)n' 6 (q; 
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(7.16) 



In this representation we have taken into account that in each chain (7.13) a fermion bubble II'-^ may be 
coupled only to a boson bubble II' 6 , while boson bubbles may be coupled to each other due to the presence 
of the superexchange term in the Hamiltonian. We have also taken into account that the terms in IT 3 ' 10 (q), 

which could give nonzero contribution to the matrix product of IT 6 ) (q) and n'^(q), vanish for q = 0. This 
allowed us to replace IT 6 ) (q) by IT ) = h/N in the leftmost and the rightmost factors of the second term 
in Eq. (7.15). 

It is convenient to use in (7.15) the basis set j^(x), n 6 '-' (x) instead of the set M^x, ■ ■ -,Mg vx , Eq. (2.7). 
Here 



(7.17) 
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while j 6 '^ (x) are given by Eq. (5.14) which we, using Eq. (2.40), rewrite as 
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With the definitions (7.17) and (7.18) our conclusions below will remain correct even when the parameters 
, 6 have nonzero imaginary parts. 
Using Eqs. (7.17), (7.18), and (2.40), the hopping term H t in Eq. (2.5) may be written as 
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H t 
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where we made the 1 — > N generalization and denoted as n v = 2l 2 t v M4 V Me v the common MF expectation 
value of the operators n 6 '-' (x) per color. Therefore, the matrix h v does not couple j b ^ (x) with n^(x), while 
the current-current matrix elements of h v are h„jbjj = h„jjjb = —l/n v . On the other hand, following the 
same route as in the derivation of Eq. (7.8), with the help of Eqs. (Al) and (A7) we obtain for the matrix 

elements of matrix 11'-^ (q) 

1 . f 

-gy{Jn(l)Jr,'(-<l))MF, 
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the following relations: 



3r,3 r i 



11%-, (q) 



-jjy (${<!) n f v ,(-q)) 



MF 



d qv iU 'Lf (q) 
iqv-n'-j-j (q) 

Jri n I 



Nn v d qv , 



(7.21) 



31 



n'^(q) is a regular function of q at q = 0, Eq. (7.21) means simply that II .^.j. (0) = Nn v S vv i and 



T'/ 



II s ~ s (0) = 0. It is clear now, that the second term in Eq. (7.15) does not contribute to the MF electromagnetic 
kernel (7.12) at q = 0, while the contribution of the first term is equal to 

e 2 n;' (0) = ?, MF <W, (7-22) 

J 1 J r,' 

where /3? IF = e 2 Nn v is the MF expectation value of p v , Eq. (7.3). Therefore, in agreement with Refs. [24, 25] 
on the MF level the system is in the normal state [26]. 

We note however, that this derivation, as well as the derivation of Ref. [25] and the functional derivation of 
Ref. [24] , is valid only on the MF level. The cancellation of the singular contributions to the MF electromagnetic 
kernel occurs between the diagrams, which differ only by one fermion RPA polarization bubble. Consider the 
higher order diagrams for the electromagnetic kernel, which have the AF propagator in the intermediate state. 
Among these diagrams one can easily find those, which can not be cancelled in this way, for example the 
diagrams which use the term ji(x), Eq. (5.16), for the current vertices in K^ e ,(q). On the other hand, there 
are no reasons for a cancellation between these diagrams, and so they produce a nonzero singular contribution 
to the electromagnetic kernel. Therefore, taking into account the corrections of higher orders in 1/N , we 
conclude that in the phase with nonzero v r , B v , F v the model exhibits superconductivity even when neither 
monons nor electrons are paired. 

VIII. Discussion 

In this paper we have presented the theory of the t-J model in the presence of the external EMF. We have 
developed the 1/N expansion for this model and shown that all the general properties of the theory are satisfied 
in every order of the 1/N expansion separately. We have also shown that superconductivity may exist in this 
model even in the absence of the electron pairing. Can our approach give a unified picture of the properties of 
the copper oxides? Let us consider the properties of different phases of the t-J model and the possible form 
of its MF phase diagram. The phases are classified first by the values of the order parameters A v , B v and F v , 
and second by whether the spinon subsystem is in the normal or superfluid state. As was discussed in Sec. 
II, the parameters B v and F v obtain nonzero values simultaneously. We will label different phases by writing 
their nonzero order parameters; the symmetric phase will be denoted as S. 

We must now distinguish between the three- and two-dimensional phases. In the three-dimensional phases 
all three components of the order parameters A v and (or) B v , F v are different from zero. In this case, if the 
spinon subsystem is in the superfluid state, it means the real Bose-condensation of spinons with nonzero order 
parameter v. We can have, therefore, the following asymmetric three-dimensional phases: A, BF , ABF , vA, 
vBF, and vABF. One can see from Eq. (2.42) for wu, that in the i>.4-phase Bose-condensation happens in 
the state with momentum h cr) = 7r/2; in the i>B_F-phase, with momentum h cr) = 0; and in the i>_4B_F-phase, 
with momentum < h cr) < ir/2 [11]. As was shown in [12], this Bose-condensation leads to the emergence 
of the long-range spin order, characterized by a presence of nonzero magnetization with the wave vector 2k c : 
(£>2k ) 7^ 0. We have, therefore, ferromagnetic long-range order (LRO) in the i>B_F-phase, Neel order in the 
i>.4-phase, and intermediate ("spiral") LRO in the i>_4B_F-phase [11]. 

In the two-dimensional phases only two components of A v and (or) B v , F v are different from zero; we 
will mark these phases with the index "2". According to Hohenberg [27], for nonzero temperature in these 
phases the spinon field cannot have nonzero expectation value; the superfluid state of the spinon subsystem 
there is the Berezinskii-Kosterlitz-Thouless state [28]. We will label these phases by the superscript a SF v for 
"superfluid". The two-dimensional asymmetric phases therefore will have such labels as ABF2, ABFrf F , etc. 

Consider now the properties of the different phases. In the S'-phase all order parameters are equal to zero. 
The gauge symmetry is not broken, and the model is in the confining state. Only gauge-invariant objects, 
like the electron, can propagate, while the propagation of the separated spinons and monons is not possible. 
On the MF level this is manifested by the momentum-independence of their Green's functions. The currents 
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j^'^(x), Eqs. (5.14) - (5.15) vanish, so on the MF level the S'-phase is insulating. There is also no magnetic 
order in it. 

In the phases A2 and A the model exhibits short- and in the phase vA, long-range magnetic order. Like in 
the S'-phase, the currents jt'^'^(x) here vanish, so on the MF level these phases are insulating. The physics of 
these phases is determined by the presence of the gapless modes - a vx and, in the i>.4-phase, spin waves. For 
the case of half-filling these phases have been treated in [12, 29], and we expect that doping does not change 
essentially their behavior. 

In the phases with nonzero B and F the currents jt'^'^(x) do not vanish. In the BF- and Bi^-phases 
the model exhibits the normal conductivity. Here the physics is determined by the gapless gauge field a vx . 
The fluctuations of the other parameters of AF 4> avx , a = 3, . . . , 6 have gaps. Therefore, their contribution to 
the AF propagator D(q) in the intermediate state in (7.11) may be neglected compared to the contribution of 
a vx , which leads to the Ioffe-Larkin composition rule [18, 24]. We also expect, that other results of [18], where 
the similar phase was treated in the slave-boson formulation, will be applicable here, which in the first place 
concerns the linear temperature dependence of the resistivity. 

The phases ABF and ABF2 are conducting with a short-range magnetic order. As was discussed in Sec. Ill, 
the gauge fields a vx and a vx here acquire gaps due to the Anderson-Higgs mechanism. The spin excitations 
here also have a gap; these are spin-gap phases. The gap in the excitation spectrum of a vx destroys the 
mechanism which produces the linear temperature dependence of the resistivity [18]. From the experimental 
point of view, the correlation between the emergence of the spin-gap and the deviation of the temperature 
dependence of the resistivity from the linear law was recently discussed in [30]. 

As was shown in the previous section, the phases vBF and vABF are superconducting. In the correspond- 
ing two-dimensional phases BFrf F and ABFrf F the vertices V v (q,k), Eq. (5.21), do not vanish. Then, 
although we have not proven it, we also expect (see Ref. [22]), that the current-current correlation function in 
these phases has a desired zero-momentum singularity. Consequently, these phases are also superconducting. 
The superconductivity exists in all these two- and three-dimensional phases without monon or electron pairing. 
We note here the essential difference between two-dimensional superconducting phases ABFrf F and BFrf F , 
and their three-dimensional counterparts vABF and vBF: in the two-dimensional case all collective modes 
have gaps, while in the three-dimensional case spin waves are gapless. This difference in the spectra must 
manifest itself in the neutron scattering. 

The possible schematic forms of the MF phase diagrams for the three- and two-dimensional scenarios 
are presented in Fig. ll(a,b). We note that the very existence of the phases shown there, and their exact 
bounds may be established, of course, only by the numerical solution of the self-consistency equations (2.37) - 
(2.39). The general structure of the phase diagram, however, is dictated by a simple condition of continuous 
dependence of the MF order parameters from the temperature and the hole density. This condition leaves a lot 
of freedom in placing the phases on the phase diagram. Therefore, a lot of different types of phase diagrams are 
possible. Schematically, the MF phase diagram for the three-dimensional scenario may look like that shown 
in Fig. 11(a). The three lines on the diagram, which may be called the i>-line, the .4-line, and the 5_F-line, 
show the boundaries of the phases with nonzero values of the corresponding order parameters. The properties 
of these phases were considered above. The different phase diagrams may be obtained by moving the lines 
on Fig. 11(a) with respect to each other. The interesting feature of this phase diagram is the existence of 
the direct antiferromagnetic-superconductor transition. Experimentally, such a transition was observed in the 
electron-doped compounds Nd2- a; Ce :c Cu04_j, and P^-^Ce^CuC^-j, [31], so we expect that these compounds 
are described by the phase diagram of Fig. 11(a). According to the discussion above, the superconducting 
phases in these materials have magnetic LRO and gapless spin excitations. It would be interesting to check 
these predictions experimentally. 

In this work we did not consider in detail the BKT-transition in the two-dimensionaH-J model. Neverthe- 
less, we expect that the schematic phase diagram for the two-dimensional scenario might look like that shown 
in Fig. 11(b). In addition to the lines presented on Fig. 11(a), we have here At and S'-F-lines. Again, by 
moving the lines on Fig. 11(b), the alternative phase diagrams may be obtained. Thus, the phase diagram of 
Fig. 11(b) contains a number of possible versions of the phase locations. Now the comparison of this phase 
diagram with the experimental phase diagrams of the hole-doped high-T c superconductors La2- a; Sr :c Cu04, 
YBa2Cu307_,5, etc. suggests that each of these superconductors may be described by one of the possibilities 
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FIG. 11. Schematic MF phase diagrams. See discussion at the end of Sec. VIII. 



provided by Fig. 11(b); we expect therefore that the hole-doped high-T c superconductors are two-dimensional 
in the sense discussed above. This conclusion agrees with the recent experimental results of Shibauchi et al. 
[32], who studied the anisotropic penetration depth in I^-^Sr^CuC^ crystals. The authors of [32] conclude 
that the electronic properties of these crystals are essentially two-dimensional. In connection with that, it 
would be interesting to obtain similar experimental data for Nd2- a; Ce :c Cu04_j, and P^-^Ce^CuO-i-j, crystals 
and to check whether they exhibit three-dimensional behavior as predicted by our approach. 

Finally, we would like to note that our approach contains in part the idea of charged bosons, which are 
responsible for the electromagnetic properties of the high-T c superconductors. The same idea is the basic 
one for the charged Bose-liquid model of the copper oxides [33]. The authors of [33] state that the unusual 
properties of the high-temperature superconductors in the superconducting as well as in the normal state 
are compatible, correspondingly, to those of a degenerate and nondegenerate state of a charged Bose-liquid. 
They claim that independently of the particular nature of charged bosons, the phenomenological concept of 
a charged Bose-liquid provides a good description of many experimentally observed properties of the copper 
oxides. This allows us to hope that our approach will be also successful in the interpretation of experiments. 
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Appendix: The derivation of (3.12) and (6.22) 

We start these derivations with the observation that the coefficients m v. satisfy the identities 
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where for a £ Mb, e a = 1 and r a = T3, while for a £ A/, t a = — 1 and r a = 1. We will also use the notations 
Th = T3, Tf = 1, and e^j = ±1. From (2.35), (2.36), and (Al) we obtain 
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where c runs over b, f and the functions F£ were defined after Eq. (5.20). In this Appendix we will make the 
trace over Nambu indices implicit; the sign of Tr will be assumed to be included in the corresponding momenta 
summation. We will also use here the formulation without the ordering field, so we have v r = in all phases. 
Now to prove (3.12) we put in Eq. (A2) u> = 0, multiply it by 



Ae^^_ f m a , v , k ^ +f 
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with a" £ He, and take a sum over k. There is no summation over rj" in Eq. (A3) and below. Shifting a 
summation variable, we obtain then for the left hand side of (A2) 

Ne'^~e c ^2 ( ? «""Vvk-q/2 ~ "Vv'k+q/2^") F k = ile cN (d qri iig a ii +d qr ,„g a „) M a »^», (A4) 

k 

where we have used (Al) and the self-consistensy equations (2.37) - (2.39); while for the right hand side of 
(A2) we obtain 
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Rewriting (A4) as 
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and comparing it with the result for the right hand side, we obtain (3.12). 

To prove (6.22), we make in (A2) the change q — > —q, multiply it by e c /2, and use the relations e a g a = g 
and e a g a = j a to present (A2) in the form 
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which is a contribution of the diagram Fig. 4(b) to Eq. (6.22). Comparing the right hand sides of (A7) and 
(6.22), we conclude that Eq. (6.22) will be proved if the contributions of the diagrams of Figs. 4(c) and 4(d) 
to Eq. (6.22) cancel each other. 

Let us show that this cancellation does really happen. The contribution of the diagram Fig. 4(d) to Eq. 
(6.22) is 

icuV^iq,^ + d. qn V< d \q,k) + d. qn V< d \q,k) 

= N E EE f ^*'-f [iujV^ b \q,k') + d. q ,V^ b \q,k') + d. qv V^ b Hq,k')] (A8) 

c'=b,f a c i k' 

xF^ + ^e- tq -^m ari , k ,W c av ,(q,k) 

(compare with the last Eq. (5.21)). Using here (A7) and shifting the summation index, we can rewrite this 
expression as 

2 N J2 J2J2 €c ' e ~ 1 ^ { m a,k'-(i/2^-T a m arikl+(l/2 ) j F^',W^(q,k) 

C '=6,/5c^*' ^ (A9) 

c'=b,f a c i k' 

where we have used (Al). Finally, using the self-consistency equations, the last expression may be written as 

" % 2 N H (d. qri ~g a + rf_ OT i«) M av W^(q, k). (A10) 

Since this is just equal to minus the contribution of the diagram Fig. 4(c) to Eq. (6.22), our proof is completed. 
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